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We predict novel phenomena in the behavior of an ultra- 
cold, trapped gas of fermionic atoms. We find that quan- 
tum statistics radically changes the collisional properties, spa- 
tial profile, and off-resonant light scattering properties of the 
atomic fermion system, and we suggest how these effects can 
be observed. 

PACS numbers: 32.80.Pj, 05.30.Fk, 32.90.+a 



The experimental ability to cool a dilute gas of atoms 
to temperatures well below one microKelvin has intro- 
duced a new area of research, highlighted by the ob- 
servation and study of Bose-Einstein condensates (BEC) 
. Interparticle interactions in these quantum degen- 
erate gases are weak and therefore amenable to theoret- 
ical treatment; in addition, the interaction strength can 
be altered experimentally [0||. A quantum degenerate 
trapped gas of fermionic atoms, i.e., atoms with an odd 
total number of constituent fermions, is relatively un- 
explored theoretically as well as experimentally. An in- 
triguing possibility is that at sufficiently low temperature 
the atoms will develop pairwise correlations, analogous to 
Cooper pairing of electrons, and undergo a phase transi- 
tion to a superfluid-like state. Theoretical studies have 
concluded that this new fermionic condensate will only 
occur at very low temperatures and will be difficult to 
observe experimentally [||. We explore a trapped gas of 
fermionic atoms that is quantum degenerate but uncon- 
densed, a regime that is interesting in its own right. 

Experimental techniques similar to those used in trap- 
ping and cooling bosonic atoms such as 87 Rb, 23 Na, or 
7 Li can be employed to trap and cool fermionic atoms, 
with possible experimental candidates including 40 K, 6 Li, 
2 H, and metastable 3 He. We focus on 40 K, which has a 
large number of trappable spin states, as an example for 
the effects discused herein. We consider 10 6 atoms in 
a cylindrically symmetric harmonic trap with a radial 
trap frequency oj r = 2?t x 400Hz and an axial frequency 
Xuj r where A = O.f. The characteristic temperature for 
quantum degeneracy, the Fermi temperature, is given by 
T F = |^(6AA0 1/3 = l.GnK f§. (Note that the criteria 
for quantum degeneracy for fermions is approximately 
the same as for bosons, since the BEC phase transition 
occurs at T c = r f±(\N/l .202) 1 / 3 ). While there is no 



phase transition for fermions at Tp, the behavior of the 
system changes below Tp where quantum statistics be- 
comes important. 

We first examine collisional properties of the trapped, 
fermion gas. In the usual BEC experiment, the atomic 
sample is cooled optically then loaded into a magnetic 
trap. To avoid losses due to spin-exchange collisions, typ- 
ically only atoms in a single spin state are loaded into the 
trap. Cooling then proceeds by forced evaporation, i.e., 
removal of the most energetic atoms and rethermaliza- 
tion of the gas via elastic collisions. For bosonic atoms 
at these low temperatures the elastic collisions are pre- 
dominately binary and s-wave in character. For fermionic 
atoms the dominant contribution must be p-wave since 
Fermi-Dirac (FD) statistics prohibits s-wave collisions 
between identical particles. However, the p-wave colli- 
sion cross-section is suppressed at low temperature due 
to the centrifugal barrier. The threshold collision en- 
ergy E t h(l) for a given partial wave I to contribute to 
the scattering can be approximated by the height of the 
centrifugal barrier, 
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Using the Ce coefficient for K-K collisions M gives a p- 
wave threshold energy Eth(l = 1) = 100/xK. Thus, below 
100 fiK, the elastic collision rate in a gas of identical 40 K 
atoms will rapidly vanish, unless there exists an acciden- 
tal p-wave resonance. 

Sympathetic cooling H can be used to circumvent 
this obstacle to evaporative cooling. If bosonic atoms 
as well as fermionic atoms are loaded into a magnetic 
trap, the bosonic atoms can be cooled evaporatively and 
the fermionic atoms will be cooled "sympathetically" via 
(s-wave) elastic collisions with the bosons. A second al- 
ternative is to load two species of fermionic atoms into 
the trap, allowing s-wave elastic collisions between non- 
identical fermions to rethermalize the gas during evapora- 
tive cooling. For the remainder of this paper, we assume 
one can successfully cool a gas of fermionic atoms to the 
quantum degenerate regime and that the background gas 
used for sympathetic cooling can be removed, without 
significantly disturbing the fermionic gas. 

The quantum statistical suppression of s-wave inter- 
actions between identical fermions makes an ultracold, 
trapped gas of fermionic atoms an excellent realization 
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of an ideal gas. This could be seen experimentally with a 
measurement of the elastic rethermalization rate jicj of a 
single component compared to a two-component (bosonic 
and fermionic atoms, or two species of fermionic atoms) 
gas. Alternatively, a study of shape excitations of the 
trapped atom cloud jllj will reveal this effect. For a 
gas of identical fermions in a harmonic trap, shape oscil- 
lations will occur at exact multiples of the characteristic 
trap frequency and will damp at a rate determined by the 
elastic collision rate in the cloud. At temperatures be- 
low the threshold for p-wave collisions (« 100/iK), which 
is still well above the temperature of quantum degener- 
acy, the ratio of the p-w&ve to s-wave elastic cross-section 
<j p /<j s scales as (ka) 4 where k is the atomic wave vector 
and a is a characteristic range of the potential. For exam- 
ple, a collision energy of 100/iK and a range a = 50 Bohr 
radii gives a p /a s « 10~ 3 . The excitation lifetimes at low 
T, therefore, will be limited by the lifetime of atoms in 
the trap (typically 100's of seconds). For comparison, 
in a two-component fermion cloud these excitations will 
have shorter lifetimes (a reasonable estimate is « 100 ms, 
the lifetime seen for shape oscillations in non-condensed, 
bosonic atom clouds |ll|]) and at sufficiently high densi- 
ties should exhibit a small frequency shift jl^] due to the 
s-wave interactions. Interestingly enough, anomalously 
low damping rates could be observed in a dilute gas of 
fermions before they are observed in atomic BEC (due 
to different physics). 

Quantum statistics also changes the collisional proper- 
ties of multi-component fermionic atom clouds. Atoms 
in two (or more) spin states could either be loaded ini- 
tially into a magnetic trap or created by an rf or mi- 
crowave pulse applied to a cold single-component cloud. 
At low temperatures (T < 100/iK) quantum statistics 
modifies the inelastic collisional loss rates by prohibiting 
s-wave spin-exchange collisions involving an initial or fi- 
nal state consisting of identical fermions. For example, 
40 K atoms in the two spin states \F = 9/2, mp = 9/2) 
and 1 9/2, 7/2} are stable against spin exchange at low 
temperatures (note the inverted ground-state structure 
shown in Fig. 1). Furthermore, at quantum degenerate 
temperatures otherwise allowed spin-exchange collisions 
can be suppressed through final state occupation. For ex- 
ample, an ensemble of K atoms in the three spin states 
|9/2,9/2), |9/2, 7/2), and [9/2,5/2) is subject to the fol- 
lowing spin-exchange collision: 9/2, 5/2) + 1 9/2, 7/2) — > 
|9/2, 3/2) + [9/2, 9/2). This process will limit the life- 
time of the 1 9/2, 5/2} atoms in the trap. However, for 
T/Tp <C 1 the Pauli exclusion principle will suppress this 
mp changing collision since the final state, 1 9/2, 9/2} with 
an appropriate energy and momentum conserving exter- 
nal state, is likely to be occupied. This suppression of 
spin-exchange could be observed as an increased lifetime 
for the 1 9/2, 5/2} atoms in the three-component cloud at 
quantum degenerate temperatures. 
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FIG. 1. Schematic of the ground-state hyperfine levels of 
40 K (shown with exaggerated Zeeman splittings). 

We now consider the interaction of atoms with light, 
which is a powerful tool for probing trapped atomic gases 
and has been studied theoretically for ultracold bosons, 
and to a lesser extent, fermions 13 - 17 . One effect of 
quantum statistics on light scattering, which has been 
essential for studying BEC with optical imaging tech- 
niques, is simply that the scattered light reflects the spa- 
tial profile of the atom cloud, which in a harmonic trap is 
directly related to the quantum statistical occupation of 
trap energy levels. For example, light scattering images 
of trapped, bosonic atoms reveal deviations from the clas- 
sical gaussian spatial profile even at temperatures slightly 
above the BEC phase transition p|p^|. The analog of 
this effect for fermions provides a straight-forward mea- 
sure of the quantum statistics, which is useful to quantify 
for comparison to other probes of quantum degeneracy. 

A trapped fermion cloud will have a larger spatial 
extent and a lower peak density relative to Maxwell- 
Boltzman (MB) particles (see Fig. 2 inset). This effect 
will be small unless T is well below T F @,||,|l[, but 
may be revealed in careful studies of the cloud profile as 
a function of temperature. For example one could ob- 
serve the momentum distribution by imaging the atoms 
after a period of free expansion from the trap jjj. The 
expected two-dimensional expanded cloud image can be 
calculated from the semiclassical momentum distribution 
for fermions in a harmonic trap |2l}| , integrated through 
in one dimension. For sufficiently long expansion times 
the initial spatial distribution can be ignored. Interac- 
tions can also be ignored for a single-component fermion 
cloud for the reasons discussed previously. A simple mea- 
sure of the extent to which FD statistics makes the ex- 
panded cloud image non-gaussian can be extracted by 
comparing a gaussian surface fit to the entire image and 
a fit to only the outer wings of the cloud image. Fitting 
the full cloud image reveals the enhanced width due to 
FD statistics (larger average energy) while a fit to the 
outer edges of the cloud gives a width that more closely 
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reflects the cloud temperature. In Fig. 2 we plot the ratio 
of the gaussian widths from these two fits, as a function 
of reduced temperature T/Tp. 
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FIG. 2. Ratio of the widths derived from gaussian surface 
fits to the full expanded fermion cloud image and to only the 
outer edges of the cloud a/a W i ngs vs. T/Tf- The wings were 
chosen to include approximately one fifth of the total number 
of atoms. The inset shows traces through the cloud image for 
FD and MB particles, calculated for 40 K, T/T F = 0.3, and a 
10 us expansion. 

A second effect of atom quantum statistics on light 
scattering comes from an enhanced or decreased proba- 
bility for scattering atoms into occupied trap levels. As- 
suming a two-level atomic system, an atom that scatters 
a photon receives recoil momentum that projects it onto 
new harmonic trap levels, which may or may not be al- 
ready occupied. (Experimentally a two-level system may 
be approximated by driving the atomic cycling transition 
with circularly polarized light whose quantization axis co- 
incides with the direction of the trap magnetic field.) For 
bosons an enhanced probability of scattering into occu- 
pied final states is predicted to give an extremely broad 
line width |l5| ]. For fermions, the Pauli exclusion prin- 
ciple does not allow scattering into occupied final states 
resulting in a blocking effect fll6| , [l7 19 which implies a 
narrowing of the linewidth. 

We explore this effect by calculating the scattering rate 
of weak, off- resonant, monochromatic light from an opti- 
cally thin cloud of fermionic atoms. This subject has also 
been treated by Javanainen and Ruostekoski jlTj who 
look at the spectrum of light scattered from a spatially 
homogeneous gas. We focus here instead on the angular 
distribution of scattered light, and treat a harmonically 
confined gas. We start with the spectral density function 
(Eq. 47 in p^] ) describing the scattered light intensity 
at r with frequency u>l + w, where ujl is the frequency of 
the probe laser, 



x j d 3 r 2 exp Mt+iAk -( ri - r2 ) 

x {^{v x , 0)V(r x , 0)^(r 2 , t)^(ra, *)), (2) 

where ?iAk is the recoil momentum of the atom, 
^-p-/(r)M(e, h) is the scattered intensity from a sin- 
gle free atom, and ip denotes the atomic ground-state 
field operator. In Eqn. 2, the excited-state field op- 
erators have been replaced by an adiabatic solution in 
terms of the ground-state field operators jl^] . We then 
use the finite temperature version of Wick's theorem p^ ] 
to evaluate the expectation value of the four field opera- 
tors, taken with respect to the atoms' many-body wave 
function. We use the unperturbed Green's function for 
fermionic atoms confined in a harmonic trap, and ignore 
coherent scattering which is relevant only for very small 
scattering angles and does not depend on the atom quan- 
tum statistics. Taking the ratio of the scattered intensity 
for FD compared to MB particles gives 



'FD 



'MB 



> w ) = ^ E Et 1 - f(E n ,T)]f(Ei,T) 



73 ^ J* 



{p + h^)^{p)\ 2 5{E n -E l -hLo). (3) 



Here i 



S(r,w) = ^J(r)^M(e,ft) 



dt / (frx 



ixiiy,iz) and n = {n x ,n y ,n z ) denote the ini- 
tial and final harmonic trap levels of the atom, f(E, t) = 
-^e- — is the FD function, and z is the fugacity. 

z exp k B T + 1 

The single particle, momentum-space wave function for 
a harmonically confined atom is denoted by </> n (p), and 
the integral above has an analytic form |23|| . Eqn. 3 can 
also be understood from Fermi's Golden Rule, where one 
integrates over all possible final states n and averages 
over initial states i. Experimentally one would detect 
the scattered photons, with the detector position defining 
the scattering angle (9, </>) relative to the incident probe 
beam (along i) and conservation of momentum deter- 
mining the atom recoil. For = 767 nm light at the 
40 K cycling transition, the recoil energy E reco . a = 

is 2\hu r . The ratio E reco u/ Ep, an important parameter 
in characterizing the light scattering properties, is then 
0.25 for our calculations. 

Figure 3 shows the angular dependence of the scat- 
tered light intensity, integrated over all scattered photon 
frequencies, for fermions relative to that for MB parti- 
cles. (For the temperatures we consider the scattering 
is essentially isotropic in </>, reflecting the isotropy of the 
momentum distribution at thermal equilibrium.) The 
blocking effect of FD statistics appears as a suppression 
of the photon scattering rate, predominantly at small an- 
gles where the small recoil momentum means that the 
final states are mostly occupied. This effect could be 
studied by monitoring the scattering rate at small angle 
as a function of temperature. Alternatively, one could 
measure the ratio of the scattered light intensity at a for- 
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ward angle to that at a backward angle, revealing the en- 
hanced backward scattering at low reduced temperatures 
p4| . In addition to modifying the angular dependence of 
the scattered light, FD statistics also suppresses the to- 
tal scattering rate (see Fig. 3). At T/T F = 0.24, for 
instance, the lifetime of the atomic excited state grows 
to 1.7 times the natural atomic transition lifetime. This 
effect is reminiscent of enhanced atomic excited state life- 
times seen in cavity experiments ]25|], although in our 
case the lifetime grows because of a suppression of atom, 
rather than photon, final states. The spectrum of scat- 
tered light, while more challenging to measure, can be 
calculated straight-forwardly from Eqn. 3. The block- 
ing effect of FD statistics shifts the spectrum to lower 
frequency as shown in the inset to Fig. 3. 
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FIG. 3. Angular dependence of scattered light intensity for 
trapped, fermionic atoms, relative to that for classical parti- 
cles, shown for E recoi i = 21h,uj r , Ef = 84hu r , and T/Tf~1.0 
(up triangles), 0.47 (circles), 0.24 (squares), and 0.10 (down 
triangles). The error on the calculated points, due to trun- 
cating the summation in Eqn. 3, is less than 0.5 %. The inset 
shows the spectrum of the scattered light for T/Tf = 0.24 
and 9 = 17°. 

The large impact of FD quantum statistics on light 
scattering can be used to probe quantum degeneracy. A 
large blocking effect requires both small T/Tf (see Fig. 
3) and E reco u less than or comparable to Ef- This second 
condition implies that the experiment is best performed 
with large N and a tightly confining trap. Another way 
to enhance the observed effect is to use a focused probe 
beam with a waist smaller than the cloud size, thus inter- 
rogating atoms in lower harmonic oscillator levels prefer- 
entially. 

In conclusion, we have examined the emergence of 
quantum statistical effects on the behavior of a trapped 
gas of fermionic atoms. Collisional properties as well as 
light scattering exhibit striking changes, such as the sup- 
pression of radiative and non-radiative decay processes, 
and provide convenient probes of quantum degeneracy. 

This work is supported by the National Institute of 
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